Preheating a bouncing universe 
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Preheating describes the stage of rapidly depositing the energy of cosmological scalar field into 
excitations of other light fields. This stage is characterized by exponential particle production due 
to the parametric resonance. We study this process in the frame of matter bounce cosmology. Our 
results show that the preheating process in bouncing cosmology is even more efficient than that 
in inflationary cosmology. In the limit of weak coupling, the period of preheating is doubled. For 
the case of normal coupling, the back-reaction of light fields can lead to thermalization before the 
bouncing point. The scenario of matter bounce curvaton could be tightly constrained due to a large 
coupling coefficient if the curvaton field is expected to preheat the universe directly. However, this 
concern can be greatly relaxed through the process of geometric preheating. 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 



Non-singular bouncing cosmologies can resolve the cos- 
mological "Big Bang" singularity of Standard Cosmol- 
ogy and have hence attracted a lot of attention in the 
literature. Non-singular bounces were proposed a long 
time ago They were studied in models motivated by 
approaches to quantum gravity such as modified grav- 
ity models 0, Q, higher derivative gravity actions (see 
e.g. [mil), non-relativist ic g ravitational actions UH, 
brane world scenarios fiol torsion gravity fl2j. "Pre- 
Big-Bang" [3 and Ekpyrotic QJ] cosmology (in which 
case the conjectured bounces are classically singular), or 
loop quantum cosmology [IH . Bouncing cosmologies can 
also be obtained using arguments from super-string the- 
ory. For example, the String Gas Cosmology scenario 
[la . may be embedded in a bouncing universe as re- 
alized in [lq . Non-singular bounces may also be stud- 
ied using effective field theory techniques by introducing 
matter fields violating certain ener gy conditions, for ex- 
ample non-conventional fluids 119 . [2C ] , quintom matter 
[2l| - [23j . or ghost condensates [2J-[26|. A non-singular 
bounce may also be achieved in a universe with open 
spatial curvature term (see e.g. [27|, H3|)- A specific re- 
alization of a quintom bounce occurs in the Lee- Wick 
cosmology studied in (29[ . Various original bounce mod- 
els were reviewed in Ref. tSfJ . 



In the context of studies of bouncing cosmologies it 
has been realized that fluctuations which are generated 
as quantum vacuum perturbations and exit the Hubble 
radius during a matter-dominated contracting phase lead 
to a scale-invariant spectrum of cosmological fluctuations 
today [UHIl (see also [33| for an earlier discussion). This 
yields an alternative to inflation for explaining the cur- 
rent observational data, which is dubbed as the matter 
bounce (see e.g. [34j for a recent review). However, in 
this scenario there are still unclear issues. For example, 
since scalar and tensor modes grow at the same rate in 



a matter-dominated phase of contraction, the tensor-to- 
scalar ratio is typically too large to be consistent with 
current observational data. Also, the possible generation 
of entropy and particles during the non-singular bounce 
phase needs to be studied in detail. In order to solve 
the first question, the scenario of matter bounce curva- 
ton was proposed in [3oT |. In this paper, we will focus 
on the question of particle production in matter bounce 
cosmology. 

The issue of particle production during the bouncing 
phase is important for various reasons. First of all, in 
the context of cold symmetric initial conditions and vac- 
uum initial conditions for the fluctuations in a matter- 
dominated phase of contraction it is important to un- 
derstand how the currently observed entropy of the uni- 
verse can arise. In this paper we will explore the possi- 
bility that this radiative entropy is generated during the 
bounce phase, in analogy to how in inflationary cosmol- 
ogy the current entropy is generated during the "reheat- 
ing" phase. Secondly, it is important to show that the 
entropy produced in the initial stages of the bouncing 
phase (i.e. before the bounce point is reached) is not 
too large to disrupt the bounce. In many models of a 
non-singular bounce (e.g. the two matter field quintom 
bounce or the Lee- Wick bounce) the bounce is unstable 
to the presence of a substantial amount of radiation [3fj| , 
since radiation red-shifts as a~ A whereas the regular and 
ghost matter fields whose interaction leads to the cosmo- 
logical bounce red-shifts at most as a -3 , where a(t) is the 
cosmological scale factor. 

We will work in the context of realizations of the 
bounce using two matter scalar fields, one with regular 
kinetic term, the other with phantom sign . A cru- 



1 For readers worried about the ghost instability of this model |37|| 
we point out that we are treating our model as an effective field 
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cial ingredient of the bounce model which we will make 
use of here is the modeling of regular matter by an os- 
cillating scalar field condensate. The oscillating scalar 
field condensate resembles the inflaton condensate dur- 
ing the preheating period during which the inflaton is 
coherently oscillating about the ground state of its po- 
tential energy. We know from studies of reheating in 
inflationary cosmology what a very interesting preheat- 
ing dynamics describes the matter transfer between the 
inflaton and regular matter. In a similar way, we expect 
that the oscillations of the scalar field matter condensate 
in the matter bounce scenario can excite fluctuations and 
thus produce radiation. This is the topic we study in this 
paper. 

Recently, the theory of preheating of the universe has 
become one of the most important issues in inflationary 
cosmology. During the preheating stage the energy is 
quickly transferred from the primordial inflaton scalar to 
matter fields which couple to the inflaton. This process 
occurs after the slow-roll period of inflation has ended, 
but earlier thermal equilibrium is established. The en- 
ergy transfer from the inflaton to matter was initially 
analyzed using first order perturbation theory and dis- 
cussed in terms of the decay inflaton quanta into Stan- 
dard Model particles [38l— [40| . However, it was realized 
that this analysis misses important effects due to the co- 
herence of the inflaton condensate [4l| . The generation 
of matter in the inflaton condensate was then studied in a 
semi-classical analysis involving non-perturbative effects 
in Ref. [4l|. It was found that a parametric resonance 
instability plays a crucial role. These parametric reso- 
nance effects after inflation were further studied in Refs. 
}42l l43j and the particle production process was then an- 
alyzed in detail in [44| . In recent years, this topic has 
been extensively studied in the literature in the frame- 
work of inflationary cosmology f45r]47j . and we refer to 
Refs. [48| and [4!| for comprehensive reviews. 

One of the interesting predictions made by preheating 
is that the metric fluctuations could be amplified due to 
entropy fluctuations [5(1 l5l|. even these fluctuations are 
outside the hubble radius. In the context of inflationary 
cosmology, it has been realized that entropy fluctuations 
can lead to an additional source of curvature fluctuations, 
a source which may in fact dominate (52| - [5r| . This is 
now know as the "curvaton mechanism" for generating 
fluctuations. Its application in matter bounce cosmology 
has been studied in Ref. (35| . 

In analogy with inflationary preheating, we expect the 
oscillations of the matter condensate to generate fluctua- 
tions in fields which matter condensate couples to. This 
is the process we study in this paper. As we show, this 



theory. There are realizations of the matter bounce which are, 
at least at the perturbative level, free from this ghost instability. 
This is the case e.g. for the cosmological bounce obtained in 
Hofava-Lifshitz gravity Q and ghost condensate bounce scenar- 
ios [MEl. 



process can indeed lead to the generation of radiation 
during the bounce. There are various constraints on this 
scenario. Firstly, the coupling of the matter condensate 
to the radiation field must be strong enough to allow for 
the preheating instability. A more stringent lower bound 
on the coupling arises if we demand that the preheating 
mechanism here has the strength to generate all of the 
radiation observed today. On the other hand, there are 
also upper bounds on the strength of the coupling. Not 
too much radiation is allowed to be generated prior to the 
bounce point - otherwise the radiation would destabilize 
the bounce. Furthermore, the back-reaction of the pro- 
duced particles is not allowed to be too strong to shut off 
the resonance via back-reaction. As we show here, there 
is a range of parameters where all of the above-mentioned 
constraints can be satisfied. 

The paper is organized as follows. In Section II we pro- 
vide the basic equations of matter bounce cosmology. In 
Section III we study the preheating process of the sim- 
plest example in the framework of matter bounce cos- 
mology. Explicitly, we study stochastic resonance with 
weak and strong coupling coefficients, respectively, and 
find the resonance is generically more efficient than that 
in inflationary cosmology. Section IV is devoted to a dis- 
cussion of the constraints on the model, and in Section 
V we show that the simplest version of preheating by 
letting the entropy field couple to the background scalar 
directly would spoil the scale-invariance of the primordial 
power spectrum in the matter bounce curvaton scenario. 
Instead, we can introduce a non-minimal coupling be- 
tween the entropy field and the Ricci scalar so that the 
curvaton field can be thermalized as well due to gravita- 
tional interactions. Section VI presents a summary and 
discussion. 

We make use of the convention m p i = 1/ yG in the 
present paper. 



II. BASIC PICTURE OF MATTER BOUNCE 
COSMOLOGY 



We will discuss matter bounce scenarios in which the 
background matter which is responsible for the matter- 
dominated phase of contraction is a scalar field <f> with 
standard kinetic term. In order to obtain a non-singular 
bounce in the context of an effective field theory in which 
space and time are described by the Einstein action, we 
need to introduce matter which violates the Null Energy 
Condition (NEC). There are various ways of achieving 
this. The details will not be important for our analysis, 
and therefore we simply use a general Lagrangian C g as 
the Lagrangian of the matter which violates the NEC. We 
assume that in the initial phases of contraction the fields 
involved in the NEC violating sector make a negligible 
contribution to the total energy-momentum tensor. In 
order to get a bouncing background cosmology, the fields 
of C g must become more important as the contraction of 
space proceeds, and eventually become equally important 
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as the standard matter at the bounce point. 

As is done in analyses of inflationary reheating, we will 
model Standard Model matter as another scalar field x 
minimally coupled to Einstein gravity R and coupled to 
4> via an interaction potential V(<f>, x)- Thus, our action 



R 



- V(ct>,x)+£ g ] ■ (1) 
As a simple example, we take the potential to be 

V{ct>, X ) = ^mV + i?W ■ (2) 



Our model is constructed such that there exists a non- 
singular bounce. The bounce point occurs at the time 
is which without loss of generality can be taken to be 
4b = 0. We choose initial conditions such that long be- 
fore the bounce point, the universe is contracting and 
dominated by a homogeneous condensate of the scalar 
field (f>. Thus, the equation of state averages to zero over 
time. During this period, the Hubble radius is decreasing 
linearly and H < 0. 

The model is constructed such that the contribution of 
C g becomes more and more important as the contraction 
proceeds. The matter component making up C g could be 
the Lee- Wick partner of cj> [2!|, it could be a condensate 
of ghost fields [26[ , or it could be terms in an effective 
action involving higher order spatial derivatives 0. In- 
dependent of the specific mechanism, at a certain point 
in the contraction phase there the Hubble radius 
will cease to decrease any further and a brief 2 bounc- 
ing phase begins. During this period H > 0. After the 
bouncing phase there comes a time when H vanishes, and 
after that the universe enters the expanding phase with 
a normal thermal history. 

Returning to the question of initial conditions, we can 
take the contracting phase to be the mirror inverse of the 
expanding phase of Standard Cosmology, i.e. containing 
both matter and radiation, with radiation dominating 
over matter for t > — t eq , where t eq is the usual time 
of equal matter and radiation for Standard Cosmology. 
If we take fluctuations to begin early in the contracting 
phase in their vacuum state, then these initial conditions 
will lead to a spectrum of curvature fluctuations at late 
times which is scale-invariant on large scales but makes a 
transition to a blue spectrum at scale which exit the Hub- 



ble radius after t 



-t ea . The observational constraints 



on this scenario were discussed in [57] . Alternatively, we 



2 The time scale of this bouncing phase (during which the equa- 
tion of state of matter crosses the phantom divide twice, first 
decreasing below w = — 1 and later - after the bounce point - 
rebounding to to > —1) is typically set by a high energy scale 
which enters into the construction of C a . 



can assume that there is no radiation at all in the con- 
tracting phase. In this case, the challenge is to produce 
the radiation which we currently observe by processes 
close to the bounce point. This is the scenario we study 
here. 

We consider the metric of a homogeneous flat 
Friedman- Robertson- Walker space-time, 



ds 2 = ~dt 2 + a 2 {t)dx 2 



(3) 



where t is the physical time, x denote the co-moving spa- 
tial coordinates and a(t) is the cosmological scale factor, 
and calculate the background evolution. The expansion 
rate of the universe, i.e., the Hubble parameter, is defined 
by H = a/a and it and its time derivative arc determined 
by the averaged Fricdmann equations, 



H z = 



H 



8irG 



(p) 



-4irG(p + p) , 



(4) 
(5) 



where p and p are the total energy density and pressure 
of the universe respectively. 
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FIG. 1: Evolution of the Hubble parameter H as a function of 
cosmic time t in a matter bounce cosmology. In the numerical 
computation, the parameters m and g in the potential were 
chosen to be m = 10 _6 m p ; and g = 10 -4 . The ghost part 
C g part taken to as that of a Lee- Wick scalar [29l ] with mass 
scale to be M = 50m. The times ti, t m , ts, and tr denote 
the moment when stochastic resonance begins, the time when 
H = 0, the bounce point, and the time of thermalization, 
respectively. The plot is in units of m v i. 

In Fig. Q]we numerically plot the evolution of the Hub- 
ble parameter in a nonsingular matter bounce cosmology 
(we use the specific model of [29| in which the Lee- Wick 
partner of (j> yields the ghost Lagrangian C g . One can 
read off that the Hubble parameter is initially negative 
which implies a contracting universe. In this phase 



(6) 



4 



where the pointed brackets indicate averaging. The back- 
ground scalar condensate <$> can be asymptotically ex- 
pressed as an oscillating function 



cj){t) ~ (^(i) sin mi , 
with a time dependent amplitude 



37rm|i| 



(7) 



(8) 



A very important lesson of this scenario is that the 
metric perturbations which originate as as quantum vac- 
uum fluctuations on sub-Hubble scales are able to cross 
the Hubble radius during the contracting phase and 
evolve on supcr-Hubble scale into classical curvature fluc- 
tuations with a scale- invariant power spectrum [32| . As 
was pointed out in [35| . isocurvature modes correspond- 
ing to massless scalar fields also obtain a scale-invariant 
power spectrum. Thus, with suitable choices on the pa- 
rameters, the matter bounce model can fit the current 
CMB observations and explain the formation of large 
scale structure of our universe [H, [58| . In addition, 
a matter bounce cosmology predicts a particular type 
of primordial non-gaussianities with a sizable amplitude 
[59( . However, without quantum excitations of light mat- 
ter fields, the universe described by a matter bounce 
would never arrive at the phase of thermal equilibrium, 
and thus there would be no radiation phase and no pro- 
cess of entropy production in this model. Note that the 
same conclusion would hold in inflationary cosmology if 
there were no couplings between the inflaton and matter 
fields. 

The analogy with inflation also tells us how to address 
this potential problem. We expect that the background 
field <f> couples to other matter fields, such as the x field 
introduced in Eq. ([T]). As is well known, by virtue of the 
coupling in Eq. @ , one obtains a resonant production of 
X particles for inflationary cosmology. One should expect 
a similar parametric resonance instability to occur in the 
matter bounce scenario. It is to a discussion of this issue 
to which we now turn. 



III. STOCHASTIC RESONANCE OF AN 
ENTROPY FIELD IN BOUNCE COSMOLOGY 

In inflationary cosmology, the theory of preheating pro- 
cess has been intensively studied in the literature. This 
process occurs ubiquitously as a period of stochastic res- 
onance in inflation models when the inflaton starts to 
oscillate about its vacuum following the phase of slow- 
roll. As we motivated above, an analogous process is 
expected to occur near the bounce phase when the uni- 
verse evolves from contraction to expansion. Thus, in the 
following we study the stochastic resonance of an entropy 
field in bounce cosmology due to the coherent dynamics 
of the 6 condensate. 



A. The Setup 

To start, we study the quantum dynamics of the en- 
tropy scalar x induced by the coherent oscillations of the 
condensate 0, first in the contracting phase, and then 
after the expanding period. Varying the action ([!) with 
respect to %, one gets the equation of motion of the in- 
homogeneous entropy field 



X + 3H X 



-x + g 2 4> 2 x 







(9) 



where we ignored the presence of metric perturbations. 
The effective mass of x is induced by its interaction with 
4> as given in the interaction potential ([2]). Note that, in 
our simplest case, the bare mass of the entropy field x is 
taken to be zero. This is a good approximation provided 
the mass of x is much less than that of the background 
field 0. 

Since the classical background solution for x evolves 
proportional to a _3//2 away from the bouncing phase (as 
shown e.g. in 35]), the investigation of parametric res- 
onance can be simplified by re-scaling the entropy field 



X 



,3/2 



X- Following the formulae developed in Ref. 
[44j , we can derive the following simpler form of the equa- 
tion of motion for the Fourier modes of the entropy field 
X k 

X k +w 2 k X k = 0, (10) 
where the frequency w\ is given by, 

The last two terms of (fTTj) comes from the gravitational 
effects which cancel each other in a matter-dominated 
phase. 

In a general case, the homogeneous solution to Eq. (fTU)) 
is composed of a negative frequency mode and a positive 
frequency mode as follows, 



X k {t) 



1 



a k (t)e- < / , »»«tt + i 8 i (t)e</ 4 »» 



dl 



,(12) 



and from (fTU|) the evolutions of these two coefficients a k 
and (3 k satisfy the equations, 



2w k 



a k ■ 



(13) 



Following the standard approach of canonical quanti- 
zation, one obtains the conjugate momentum Tl k for the 
field variable X k : 



n. 




(14) 
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where we have used the equations (|13[) . Imposing the 
canonical normalized commutation relations 



[X fe ,n*/] = iS(w k - Wy 



(15) 



one finds that the coefficients and j3 k coincide with 
the coefficients of the Bogoliubov transformation of the 
creation and annihilation operators in curved space-time. 
Thus, we impose the normalized initial conditions, 



i , p k = o , 



(16) 



at the beginning of the preheating period U. Moreover, 
the occupation number density for \ particles with co- 
moving wave number k is given by 



2 



in* 



\x k 



(17) 



where the last term i on the right hand side is a c-number 
used to cancel the zero-point energy. 

There is a parametric resonance instability in the above 
cosmological system which can be described in terms 
of successive periodic scatterings of the Yk modes in a 
parabolic potential, as first studied in j44|. The reso- 
nance is effective if the condition 



\gm P i\ > y to 



(18) 



is satisfied (which comes from demanding that q > 1 is 
realized at some point during the oscillatory phase of <f> 
- see the discussion below). In this case the total num- 
ber density of the x particles created during stochastic 
resonance can be estimated to be 



n x (t) 



2ira 



n k {t) 

£,3g2mr ; (t-t;) 



647r 2 a 3 ^Jirmrjit — ti 



2tt 2 



(19) 



where ti denotes the moment when stochastic resonance 
begins (this time is shown in Fig. [I]). 

The value of ti in bounce cosmology is different from 
that in inflation. In inflationary cosmology, the preheat- 
ing phase starts which <j) takes on the largest value it 
has during its phase of oscillation, whereas in the matter 
bounce the preheating phase starts at the lowest value 
of <fr when the efficiency condition for resonance is satis- 
fied (this is the condition q > 1 to be discussed below). 
Explicitly, ti = for inflation while 



U = - 



grripi 



(20) 



in bounce cosmology as will be explained in the following 
section. 

Resonance takes place for all values of k in certain res- 
onance band. In the case of stochastic resonance (which 
is realized here) , the dominant resonance band is the long 



wavelength band whose half width fcj at the initial time 
ti is given by 



k 

a i 



gm\4>i 



(21) 



Its physical wave number approximately equals to 
*/ gmrripi/5 in inflationary cosmology but m in bounce 
cosmology. The parameter rj is the largest value of the 
Floquet (growth) index within the instability band, and 
its effective value is 

r] ~ 0.18 - 0.1. (22) 



B. Broad Resonance with Weak Coupling 

Recall that the cosmological scalar (f> oscillates around 
its vacuum and the universe evolves as matter dominated 
one in the contracting phase. Its background evolution 
is described by Eqs. (j6|), (J7|) and the amplitude is given 
by ©• It is useful to use, instead of cosmic time t, the 
number of oscillations of the homogeneous scalar 0, which 
is (beginning the count at the time ti when resonance first 
becomes effective) 



N(t) 



m(t - U) 
2^ 



(23) 



In the analytic computation of a process of stochastic 
resonance, it is convenient to define a resonance param- 
eter 



g 2 J> 2 



m- 



(24) 



Studies of reheating in inflationary cosmology [44|, [49[ 
have taught us that in a dynamical space-time only the 
broad resonance regime with q > 1 can lead to a signif- 
icant enhancement of the generation of the \ particles 
and of the corresponding fluctuations on super-Hubble 
scales. Since in the framework of bounce cosmology the 
universe starts its evolution from a very low energy state 
in a contracting phase, the value of q is initially much 
smaller than 1 which indicates broad resonance cannot 
take place in most of this regime. However, along with 
the contraction, the amplitude of the background scalar 4> 
increases gradually, and consequently, q will approach to 
the boundary q = 1 and then cross over. That moment 
corresponds to the starting time ti of broad resonance 
in bounce cosmology, which - if wc follow the solution 
obtained in Eq. ([8] - is given by 



U 



gm.pi 



3tttti 2 



(25) 



with & = 2). 



As long as we can neglect the back-reaction of the 
X particles produced during the resonance on the back- 
ground tfi condensate, the condensate will oscillate with 



6 



increasing amplitude <fr until it reaches the value 
~ _ 2m p i 



which occurs at the time 



2m 



(26) 



(27) 



in the contracting phase, and then the universe will enter 
the bouncing phase with the broad resonance being re- 
placed by a process of gravitational resonance as studied 
in Sec. V. After the bounce, there is one another period 
of broad resonance in the expanding phase, which is a 
reversal of the process in the contracting phase. 

Therefore, during the process of broad resonance with 
weak coupling in bounce cosmology, the scalar <f> experi- 
ences the following number of oscillations 



N 



gm p i 
3irm 



(28) 



which is double of that obtained in usual inflationary 
cosmology. For example, for the set of parameters m = 
10~ 6 m p i and g = 10~ 4 , the above estimate gives N ~ 10. 
Correspondingly, the increase in the occupation number 
density for the k = 0.1m mode of x is given by 



A In n fe ~ 47T77fcAr-ln2 ~ 15 
if we use the value r\u ~ 0.13. 



(29) 




FIG. 3: Evolution of the variable Xk during the process of 
broad resonance with weak coupling. The horizontal axis is 
cosmic time. The initial values of the background parameters 
are the same as in Figure [2] 



results in Figs. [2] and |3l We find that the analytic es- 
timate is in good agreement with the numerical results 
shown in these figures. One can conclude that the main 
mechanism of broad resonance with weak coupling is sim- 
ilar to that in inflationary cosmology, but that the oscil- 
lation number is doubled since there are now two time 
periods of resonance as a consequence of the background 
evolution. 
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FIG. 2: Evolution of the occupation number density nk dur- 
ing the process of broad resonance with weak coupling as a 
function of cosmic time (horizontal axis). The time axis is 
displayed in units of m~ . The parameters m and g were 
chosen to be m = 10~ 6 m p ( and g = 1CP 4 . In the compu- 
tation, the mode was chosen as k — 10~ 7 m p i with resulting 
initial conditions Xk — 0.7386 and Xk = 0. 

To get a better understanding of this process, we solve 
the perturbation equation (|12p numerically and show the 



C. Pre-Bounce Resonance with Strong Coupling 

As we have analyzed in the above subsection, when 
g > \f2nm/m p i, our model will experience a period of 
broad parametric resonance as occurs in inflationary pre- 
heating. Note that we have so far neglected the back- 
reaction due to the presence of x particles. Once the 
value of the coupling constant g is sufficiently large, this 
back-reaction can no longer be neglected. This back- 
reaction can change the structure of the resonance, and 
it can also change the background evolution. We will 
study the first issue during the pre-bounce resonance in 
the following. 

The rough criterion for the back-reaction to have a 
negligible effect on the structure of the resonance is to 
demand that the change AmJ in the mass of as a con- 
sequence of interactions with the produced x particles be 
smaller than its intrinsic mass square m 2 . In the Hartree 
approximation, the change in the mass of the background 
scalar <p due to the back-reaction of x quanta is expressed 
as, 

Ami = 9 2 (X 2 ) ■ (30) 
Here, (x 2 ) indicates the quantum expectation value of 
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the square of the entropy field, which in turn is given by 



r 



dk- 



27r 2 a 3 



(31) 



The preheating process in the contracting phase can be 
divided into two periods. In the first the level of excita- 
tion of x particles is low enough to be able to neglect the 
back-reaction effects, and the frequency of oscillations of 
the scalar <f> is still determined by the mass m. Recall 
that we used the occupation number density (fT9|) to de- 
scribe the production of \- As can easily be seen, the 
expectation value pip is related to the number density 
of x particles via 



(x 2 



(32) 



where fej is the half width of the dominant resonance band 
and which before the bounce corresponds to a physical 
wave number (k/a)i ~ m. As the generation of x parti- 
cles, the back-reaction effect starts to be important, lead- 
ing to an end of the first stage of resonance at a moment 
t p i when 



n x (t pl ) 



which yields 



t 



pi 



my) 

2rjm 



If \ 
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9 2 V 2 m^ 2 9 r,m pl 



327r 3 m2 



(33) 



(34) 



where W is a Lambert W-function. 

Consequently, the oscillation number of <p in the first 
stage of preheating can be estimated by the following 
approximate relation: 



pi 



gmpi _ W{y) 
6irm 47T7/ 



(35) 



In the first stage of pre-bounce resonance, the amplitude 
of the scalar <fi increases and thus leads to q 1. 

The preheating process will enter a second stage if 
t p i < t m , i.e. if back-reaction becomes important be- 
fore the amplitude of <fi has reached its maximal value. 
This gives a criterion for the coupling constant g: a sec- 
ond stage of preheating will occur if g is larger than the 
following critical value: 



15m 47?»m^ , 
9c ^ W( T^e' 



') 



(36) 



There is one further condition on g: if g is too large, 
then back-reaction cannot be neglected at the initial time 
when the resonance condition becomes satisfied. This 
will occur if t p i > ti (recall that both of these times are 
negative) . We will call the related critical value of g by 
gi. Thus, we require g < gi with 



47T 5 



\VmpiJ 



(37) 




9 



FIG. 4: The number of oscillations N p i in the first stage of 
pre-bounce resonance as a function of the coupling constant 
g. In the calculation we took m = 10 -6 m p i and rj = 0.13. 



If we choose m = 10 _6 m p ; and r\ = 0.13, we find that 
in order to have successful preheating beginning in a first 
stage during which the x back-reaction is negligible and 
leading on to a second stage when the x particles domi- 
nate the mass of 0, we require the coupling parameter g 
to be bigger than g c ~ 2.5 x 1CP 4 but less than gi ~ 1.5. 

We have performed a numerical estimate of the oscil- 
lation number N as a function of g and the results are 
displayed in Fig. |4j One may notice that there are no 
oscillations of the scalar <f> when g is larger than of order 
0(1). This is understandable since g > 1 corresponds to 
a very strongly coupled system and so the present pertur- 
bative treatment is no longer valid. We will not consider 
such large values of g any further in this work. 

Finally, let us estimate the amount of preheating in 
the second stage, during which the frequency of the <fi 
oscillations is no longer determined by m but instead by 
5Xrms- Therefore, the resonance parameter is effectively 
modified to be 



q e ff 



2 

rms 1 



(38) 



The stage terminates at the moment t P 2 when the ef- 
fective resonance parameter decreases to 1. As a con- 
sequence, similar to what is done in the analysis in in- 
flationary preheating, we obtain the following oscillation 
number in the second stage of pre-bounce resonance: 



N, 



l>2 



ln(2V2g r \) 



(39) 



and correspondingly the number density of x particles is 
given by 



Xp2 



(40) 



Inserting the values of parameters chosen in Fig. [4j one 
finds that N P 2 < 1 in most of the allowed regime for 
g. This result indicates that once the back-reaction of x 
particles dominates over the background universe in the 
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contracting phase, there is no longer a broad resonance 
band after the bounce. 

To briefly summarize: in this section we have analyzed 
the parametric resonance of the entropy field \ arising 
in the simplest model in the framework of a nonsingular 
bounce cosmology. We find that there exist two differ- 
ent resonance regions (in time): one is a broad resonance 
region with weak coupling which starts in the contract- 
ing phase once the amplitude of the oscillating matter 
condensate field <f> exceeds some critical value, and fi- 
nally ends in the expanding phase once the amplitude 
falls again below this critical value. The total oscillation 
number during this broad resonance interval is twice that 
obtained in inflationary cosmology. Around the bounce 
point there may be a second region of broad resonance 
with a relatively strong coupling. This second region 
arises if the back-reaction of the entropy particles on the 
background condensate becomes important. If this back- 
reaction becomes important, then the existence of the 
part of the first resonance interval after the bounce will 
be affected. 



IV. CONSTRAINTS ON g 

There are two additional constraints on the model pa- 
rameters which we have not yet discussed. First, the 
resonance has to be weak enough such that the \ parti- 
cles (whose energy density grows as a -4 and thus behave 
as radiation) do not dominate the energy density at the 
bounce time (bounce time computed without taking into 
account the back-reaction of the x particles) . If they were 
to dominate, the bouncing background would be desta- 
bilize and the universe would collapse to a Big Crunch 
singularity instead of bouncing. Secondly, if we want 
preheating to explain the origin of the current entropy, 
then we require that at the end of the resonance period 
in the expanding phase the energy density in \ particles 
becomes dominant. In this section we discuss these and 
the other conditions on the coupling constant g. 

First of all, we recall from (fT5)) that no broad resonance 
at weak coupling can take place unless 



of mi we obta 



.9 > 



37TTO 

2m p i 



(41) 



If we want to back-reaction of \ particles not to cut off 
the resonance before the bounce point, we have an upper 
bound on g 



9 < 



15m 
Arjm p i 



4r/5ml, 2 
W r^e = 



7T77 



(42) 



Next, we turn to the energy density conditions men- 
tioned at the beginning of this section. The energy den- 
sity of the entropy field x can easily be obtained by in- 
tegrating the number density times mode energy of x 
modes over the dominant instability band. Making use 



2mrj(t-ti) 

Py(t) = pAti)( — ] ) — , 



(43) 



where a% denotes the scale factor of the universe at the 
beginning moment of preheating. Note that, the second 
term of rhs of (|4"3")) is arisen from the redshift effect, and 
the third term is due to the parametric resonance effect 
of preheating. 

Now we compare the energy density of entropy field 
and that of background scalar at the moment t m when 
the background scalar ceases its oscillations in the con- 
tracting phase. We require p x < p^ at this moment, and 
thus obtain an upper bound on the coupling parameter, 
which takes the form 



.9 < 



13m w ( ^ 16/13 



12/13 

V • 



Ar/nipi 



12/13 



(44) 



Next, we calculate p x and p^, at the end of the para- 
metric resonance time interval (time t/) in the expand- 
ing phase (neglecting the possible back-reaction of the x 
particles on the background). If p x > p^ at tf then the 
universe will enter a radiation-dominated period after the 
bounce. As a consequence, we get a lower bound on the 
coupling parameter, which is 



.9 > 



3m W( l ^ /2 r)m p i 



irjmpi 



3m 



(45) 




1.0x10 2.0x10 3.0x 

m 

FIG. 5: The constraints on the parameter space of the cou- 
pling constant g. In the calculation we took g — 0.13 and the 
mass of background field m is in units of m p i . 

Combining (gTJ), (|g2"|). ([H]) and (gSJ), we find that the 
parameter space of the coupling constant g is tightly con- 
strained by these inequalities. Specifically, in Fig. [5] we 
plot the contours constraining the viable parameter space 
for g In the region above the red line, g is too strong 
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so that the backreaction of \ dominates over the back- 
ground before the bounce point. The requirement of con- 
trollable growth of x particles before the bounce imposes 
an upper bound on g which is shown by the green line as 
analyzed in (|44|) . The blue line corresponds to the lower 
bound on g obtained by demanding that the universe be 
dominated by radiation after the preheating phase. Fi- 
nally, the regime below the purple line does not allow the 
occurrence of parametric resonance. 

V. GEOMETRIC PREHEATING AND THE 
SCENARIO OF MATTER BOUNCE CURVATON 

The same two scalar field model as studied here was 
used in Ref. HU to address one of the potential prob- 
lems of the matter bounce scenario, namely the fact that 
a simple matter bounce typically produces a ratio of ten- 
sor to scalar metric fluctuations which is of order unity 
and already in marginal conflict with data. The idea is to 
introduce a light field x to seed isocurvaturc fluctuations. 
The infrared modes of the isocurvaturc fluctuations can 
obtain a kinetic amplification and are converted to curva- 
ture perturbations in the bouncing phase. This so-called 
matter bounce curvaton scenario requires that the coeffi- 
cient g of the coupling between the curvaton field and the 
background scalar has to be very small to ensure that the 
isocurvature fluctuations obtain a nearly scale-invariant. 
In fact, the condition on g derived in [35| is the opposite 
of the criterium (|T8l) derived here which is required to 
obtain broad resonance. 

Fortunately, these two different pictures can be unified 
by introducing a non-minimal coupling for the curvaton 
field. From the point of view of quantum field theories 
in curved space-time, non-minimal coupling terms often 
arise to cancel the divergencies which arc widely encoun- 
tered in the ultraviolet limit of renormalization group 
flows of cosmological scalars. As a consequence, we ex- 
pect an existence of a non-minimal coupling for the curva- 
ton field in the matter bounce can provide an alternative 
channel for parametric resonance. The same idea applied 
to inflationary cosmology has been analyzed in |60l . Itnj 
and was called "geometric preheating" . In the following 
we explore the possibility of geometric preheating in the 
matter bounce curvaton scenario. 

We consider the curvaton field x non-minimally cou- 
pling to the Ricci scalar R through the following term, 

Vr = \^ 2 RX 2 ■ (46) 

In the matter-dominated contracting universe we have 

R = 6H+12H 2 

c 8^^, (47) 
m P l 



the frequency of the entropy field in Fourier space is given 
by 

wl * ^ +ff V + 8V^ 2 - (48) 
a m P i 

One notices that, even the coefficient of direct coupling 
between the curvaton x an d the background scalar <f> is 
negligible, parametric resonance is still able to be taken 
place if the effective resonance parameter is larger than 
unit at the moment of U as studied in the previous sec- 
tion. The condition for geometric preheating to arise in 
the matter bounce curvaton scenario is 

(49) 

Comparing this with the criterion of geometric preheat- 
ing in inflation (which is n > ~ I), we see that the 
requirement in our model is much relaxed. 

Note that the above estimation is obtained based on 
the following approximation: we only considered the 
time-averaged effect of the background evolution, and ig- 
nored the back-reaction of both the entropy field and the 
metric fluctuations. Thus, our result only give a quali- 
tative guideline. It is necessary to do a detailed quan- 
titative analysis of this scenario. Moreover, a general 
model involving non-minimal coupling to background 
fields might be interesting 3 . We would like to leave these 
issues to followup work. 



VI. CONCLUSIONS AND DISCUSSION 

In this paper we have studied the theory of preheating 
in the framework of a nonsingular bouncing cosmology 
with a matter-dominated contracting phase. Important 
in our analysis is that the contracting phase is dominated 
by an oscillating scalar mater field condensate. We have 
seen that preheating is driven by a period of stochastic 
parametric resonance and can in principle explain the 
origin of the post-bounce entropy. Our analysis does not 
depend much on the precise mechanism which realizes the 
non-singular bounce since the preheating takes place be- 
fore the terms yielding the non-singular bounce become 
very important. 

The epoch of preheating takes place near the bounc- 
ing phase after the amplitude of the condensate has had 
time to grow to a sufficiently large value for the paramet- 
ric resonance parameter q to become larger than q = 1. 
Depending on the magnitude of the coupling coefficient 
between the condensate and the entropy field x, preheat- 
ing will proceed either throughout the whole pre- and 
post-bounce intervals during which q > 1, or it will end 



where we have made use of time averaging and ignored 
the back-reaction of the entropy field. As a consequence, 



3 Non- minimal matter bounce cosmology was studied in |62| . 
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before the bounce. When the coupling is weak, the pro- 
cess of stochastic resonance is similar to that in inflation, 
but the oscillation number is doubled. Therefore, com- 
pared with inflationary cosmology, there will be more 
particles of the entropy field excited in a bouncing cos- 
mology. In the case of strong coupling, the back-reaction 
of the entropy field cannot be neglected and leads to an 
earlier cutoff for the time interval of stochastic resonance. 
For a typical parameter values of our model, it is found 
that parametric resonance terminates before the universe 
arrives at the bounce point. 

We also have studied an extended picture of preheat- 
ing in which the entropy field couples to the background 
scalar field mainly gravitationally. Explicitly, the entropy 
field couples to the Ricci scalar non-minimally and thus 
can also give rise to parametric resonance. In this ap- 
proach, we are able to accommodate the entropy gen- 
eration studied here with the matter bounce curvaton 
scenario [Hj]. According to our qualitative estimates, the 
condition for the entropy field to preheat geometrically 
in a matter bounce scenario is slightly easier than in the 
case of inflation. 

Note that we have not addressed the interesting ques- 
tion of how thermalization after preheating takes place 
in the context of a nonsingular bouncing universe. This 
thermalization process involves a lot of non-perturbative 
effects, and thus calls for corresponding non-perturbative 
analyses. Some numerical studies in the framework of 
inflationary cosmology were performed in recent years 
[63| - |65| . and it was found that the evolution of the en- 
tropy field enters a regime of turbulent scaling [6^. We 
expect that a similar phase could also occur in the matter 
bounce model. However, the most distinctive difference 



between preheating in inflation and matter bounce is the 
initial state. We expect that this could yield new results 
for thermalization of a bouncing universe. 

To conclude, we would like to highlight the importance 
of our analysis. The study of preheating has a lot of 
applications to other topics, e.g. to topological defect 
production [67], [Hf, primordial magnetic fields (69l - [7l| . 
induced non-Gaussianities [73 - l76| . preheating with non- 
standard kinetic terms ffl\ . and so on. The mechanism 
of preheating in bounce cosmology, since the initial con- 
dition is modified when compared to inflationary pre- 
heating, could provide a new window to explore early 
universe phenomenology combined with particle physics. 
This process is rather robust and ought to be considered 
in all bounce models. 
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